Screening effects on the excitonic instability in graphene 
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We investigate the excitonic instability in the theory of Dirac fermions in graphene with long- 
range Coulomb interaction. We analyze the electron-hole vertex relevant for exciton condensation 
in the ladder approximation, showing that it blows up at a critical value of the interaction strength 
a = e 2 /A-KVF sensitive to further many-body corrections. Under static screening of the interaction, 
we find that taking into account electron self-energy corrections increases the critical coupling to 
q c ~ 2.09, for a number N = 4 of two-component Dirac fermions. We show that the dynamical 
screening of the interaction has however the opposite effect of enhancing the instability, which turns 
out to develop then at a c « 0.99 for N = 4, bringing the question of whether that critical value can 
be reached by the effective coupling in real graphene samples at the low-energy scales of the exciton 
condensation. 
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Introduction. — The discovery of graphene, the mate- 
rial made of a one-atom-thick carbon layer, has attracted 
a lot attention as it provides the realization of a system 
where the electrons have conical valence and conduction 
bands, therefore behaving at low energies as massless 
Dirac fermions (JQ- This offers the possibility of em- 
ploying the new material as a test ground of fundamen- 
tal concepts in theoretical physics, since the interacting 
electron system represents a variant of strongly coupled 
quantum electrodynamics (QED) affording quite unusual 
effects 0-0| ■ 

A remarkable feature of such a theory is that a suf- 
ficiently strong Coulomb interaction may open a gap in 
the electronic spectrum. This effect, already known from 
the study of QED [8(, was first addressed in graphene 
in the context of the theory with a large number N of 
fermion flavors [9T-fl2|. The existence of a critical point 
for exciton condensation was also suggested from second- 
order calculations of electron self-energy corrections [13|. 
More recently, Monte Carlo simulations of the field the- 
ory have been carried out on the graphene lattice [3 EH; 
showing that the chiral symmetry of the massless theory 
can be broken at the physical value N = 4 above a crit- 
ical interaction strength a c « 1.08 [l4| . The possibility 
of exciton condensation has been also studied in the lad- 
der approximation [l6l - [T9| . leading in the case of static 
screening of the interaction to an estimate of the critical 
coupling a c ~ 1.62 for N = 4 [l6j]. Lately, the resolution 
of the Schwinger- Dyson formulation of the gap equation 
has revealed that the effect of the dynamical polarization 
can significantly lower the critical coupling for exciton 
condensation, down to a value a c w 0.92 for N = 4 20]. 

In this paper we take advantage of the renormaliza- 
tion properties of the Dirac theory in order to assess 
the effect of different many-body corrections to the ex- 
citonic instability. In this respect, the renormalization 
of the quasiparticle properties can have a significant im- 
pact, mainly through the increase of the Fermi velocity 
at low energies pli [22| . Thus, we will consider the renor- 
malization of the electron-hole vertex accounting for the 
exciton condensation in the ladder approximation, sup- 
plemented by self-energy corrections to electron and hole 



states. This dressing of the bare quasiparticles will have 
the result of increasing the critical coupling at which the 
excitonic instability takes place, going in the same direc- 
tion as the effect of screening the Coulomb interaction. 
We will see however that, incorporating the dynamical 
polarization in the ladder approximation, the screening 
effects are softened at N = 4, leading to values of the crit- 
ical coupling below those corresponding to the nominal 
interaction strength in isolated free-standing graphene. 

We consider the field theory for Dirac quasiparticles 
in graphene interacting through the long-range Coulomb 
potential, with a Hamiltonian given by 



H = iv F I d 2 r 4>i{r)j • V^i(r) 

+ £/ dVl / d2r2P(ri) F^ 
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■P(r 2 ) 
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where {ipi} is a collection of N/2 four-component Dirac 
spinors, ip i = ipjjo, and p(r) = ^(r)7o^»(r). The ma- 
trices 7 CT satisfy {7^,7^} = 2 diag(l, — 1, —1) and can be 
conveniently represented in terms of Pauli matrices as 
70,1,2 = (03, C3CT1, ozori) ® 03, where the first factor acts 
on the two sublattice components of the graphene lattice. 

Our main interest is to study the behavior of the vertex 
for the operator p m (r) = tp(r)ip(r). This gives the order 
parameter for the exciton condensation, and the signal 
that it gets a nonvanishing expectation value can be ob- 
tained from the divergence of (Tp m (q, t)p m (— q, 0)}. The 
singular behavior of this susceptibility can be traced back 
to the divergence of the vertex for (p TO (q)^(k+q)^(k)) 
at q — >• 0. We will denote the vertex in this limit (setting 
also the corresponding frequency co q — 0) by r(k, ojf.). In 
the ladder approximation, depicted diagrammatically in 
Fig. [TJ the vertex is bound to satisfy the equation 



r(k, uk) = 70+ 



d 2 p du)p r(p,Wp) 



2tt v 2 F \) 2 



V(k~p,iuj k - 



(2) 

where V(p,ui p ) stands for the Coulomb interaction. We 
will deal in general with the RPA to screen the potential, 
so that V(p,uj p ) — e 2 /(2|p| + e 2 x(p, w p )), in terms of the 
polarization \ f° r N two-component Dirac fermions. 
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The expressions (|4| and §Q are singular in the limit 
e — > 0. The most convenient way to show that all the 
divergences can be renormalized away is to resort at this 
point to a perturbative computation of T(k, 0). 
The solution of ((4]) can be obtained in the form 



FIG. 1: Self-consistent diagrammatic equation for the vertex 
of (p m (q)ip(k + q)V' t (k)), equivalent to the sum of ladder 
diagrams built from the iteration of the Coulomb interaction 
(wavy line) between electron and hole states (arrow lines). 



Eq. ([2J is formally invariant under a dilatation of fre- 
quencies and momenta, which shows that the scale of 
r(k, ujk) is dictated by the high-energy cutoff A needed 
to regularize the integrals. The vertex acquires in general 
an anomalous dimension 7,;, 2 , which governs the behavior 
under changes in the cutoff[23| 



r(k,Wfc) ~ A 7 ^ 2 



(3) 



We recall below how to compute 7^2 , showing that it di- 
verges at a critical value of the interaction strength. This 
translates into a divergence of the own susceptibility at 
momentum transfer q — > 0, providing then the signature 
of the condensation of p m (r) = ip(r)ip(r) and the conse- 
quent development of the excitonic gap. 

Self-energy corrections to ladder approximation. — We 
deal first with the approach in which electrons and holes 
are dressed by self-energy corrections, while the Coulomb 
interaction in © is screened by means of the static RPA 
with x(PjO) = (N/16)\p\/vf ■ It is known that graphene 
remains a conventional Fermi liquid even at the charge 
neutrality point, with a quasiparticle weight that does 
not vanish at the Fermi level[24|]. The most important 
self-energy effect comes from the renormalization of the 
Fermi velocity at low energies [25|, and this is the feature 
that we want to incorporate in our analysis, identifying 
vf in Eq. (f5|) with the Fermi velocity dressed by self- 
energy corrections. 

The electron self-energy corrections, as well as the 
terms of the ladder series, are given by logarithmically 
divergent integrals that need to be cut off at a given scale 
A. Alternatively, one can also define the theory at spatial 
dimension d — 2 — e, what automatically regularizes all 
the momentum integrals. Eq. ([2]) then becomes 
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where Vp(p) is the Fermi velocity dressed with the self- 
energy corrections, k = 1 + Ne 2 /32vF, and e 2 , is related 
to e 2 through an auxiliary momentum scale p such that 

p 1 - n e p 2 

e — p e . 

In the ladder approximation, the Fermi velocity gets 
a divergent correction only from the "rainbow" self- 
energy diagram with exchange of a single screened 
interaction [25J. The dressed Fermi velocity becomes 



v F (p) = v F 



1 
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with Ao = e\j 'Attkvf ■ Each term in the sum can be ob- 
tained from the previous one by noticing that 
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At each given perturbative level, the vertex displays then 
a number of poles in the variable e. The crucial point 
is that these divergences can be reabsorbed by passing 
to physical quantities such that vp = Z v (vF) Ten and 
ipip — Z^[il)ijj) Ten (the scale of the single Dirac field is 
not renormalized in the ladder approximation). 



The renormalized vertex r re 



Z^,iY can be actu- 



ally made finite with a suitable choice of momentum- 
independent factors Z v and Z^i . Z v must be chosen to 
cancel the 1/e pole in ([5]), and it has therefore the simple 

On 



structure Z v = 1 + bi/e, with b\ 



~ /16ttk(vf)t 



the other hand, we have the general structure 
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The position of the different poles must be chosen to 
enforce the finiteness of r rcn = Z^iY in the limit e — > 0. 
The computation of the first orders of the expansion gives 
for instance the result 



ci(A) 



1 



A - i log(2) A 2 - (k 2 + 120 log 2 (2)) A 3 



IOtt 2 log(2)+688 log J (2) + 15C(3) >4 
6144 A 



137r 4 +20647T 2 log 2 (2) + 144(7161og 4 (2)+371og(2)C(3)) 5 



737280 



16" ■ M ~ B V-,'. ■ 18432 (5^ + 744 l0 S 2 ( 2 )) A ' 

110^ 2 log(2)+8592 log 3 (2) + 135<(3) x 5 
+ 184320 



ca (A) = ^A 2 + ilog(2)A 3 



A 5 + 



c 3 (A) 
c 4 (A) 



= -76^8(2) A 4 
log(2) A 5 



184320 



360 log 2 (2)) A 5 



7680 



(9) 



where the series are written in terms of the renormalized 
coupling A = p~ £ Z v \o 

The physical observable in which we are interested is 
the anomalous dimension 7^,2 . The change in the dimen- 
sion of F ren comes from the dependence of Z^ on the 
only dimensionful scale p, being 7^,2 = (p/Z^dZ^/dp 
[23l | . The bare theory at d ^ 2 does not know about the 
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arbitrary scale p, and the independence of Ao = p e X/Z v 
on that variable leads to 



dX 
' 'dp 



-eX - Xbx(X) 



The anomalous dimension is thcn|2(: 



dX d log Z^ 
% dX~ 



= -Xc' 1 (X) 



(10) 



(11) 



In the derivation of (|11[) . it has been implicitly assumed 
that poles in the variable e cannot appear at the right- 
hand-side of the equation. For this to be true, the 
set of equations c' n+1 = c n c! x — b\c' n must be satisfied 
identically [26|. Quite remarkably, we have verified that 
this is the case, up to the order A 7 we have computed 
the coefficients in ©. This is the proof of the renormal- 
izability of the theory, which guarantees that physical 
quantities like 7^2 remain finite in the limit e — > 0. 

From the practical point of view, the important result 
is the evidence that the perturbative expansion of ci(A) 
approaches a geometric series in the A variable. It can 
be checked that the coefficients in the expansion grow 
exponentially with the order n, in such a way that 



ci(A) >^2d n X n 

71=1 



(12) 



A lower bound for d can be obtained from the first orders 
in ci(A). This estimate tends to increase as it is made 
from higher orders in the expansion. The assumption of 
scaling with the order n allows us to estimate a radius 
of convergence A c « 0.49 (to be compared with the value 
in the approach neglecting self-energy corrections, which 
leads to A c m 0.45 [Hj], in close agreement with the re- 
sult of Ref. [HI). The critical coupling in the variable 
A = a/n can be used to draw the boundary for exciton 
condensation in the (N, a) phase diagram, represented 
in Fig. For N = 4, we get in particular the critical 
coupling a c w 2.09, significantly above the critical value 
that would be obtained from the radius of convergence 
without self-energy corrections (a c ~ 1.53). 

Dynamical screening in the ladder approximation. — 
In the framework of the ladder approximation, one can 
also study the effect of the dynamical screening of the 
Coulomb interaction. We can go beyond the static RPA 
by taking into account the full effect of the frequency- 
dependent polarization, which for Dirac fermions takes 

the form x(p,w p ) = (7V/16)p 2 / ^Jv 2 F p 2 - w 1 [Hj]. This 

expression can be introduced in Eq. @ to look again for 
self-consistent solutions for the vertex T(k, w^). While 
this problem does not afford an analytic approach of the 
type shown before, one can resort to numerical methods 
for the resolution of the integral equation. In this pro- 
cedure, we find again that there is a critical coupling at 
which r(k, ujk) blows up, marking the boundary between 
two different regimes where the integral equation has re- 
spectively positive and negative (unphysical) solutions. 
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FIG. 2: Phase diagram showing the boundary of the phase 
with exciton condensation (m 7^ 0), obtained in the ladder ap- 
proximation with and without electron self-energy corrections 
(upper and lower dashed lines) and with dynamical screening 
of the Coulomb interaction (full line). 



In practice, we have solved the integral equation ([2]) 
by defining the vertex in a discrete set of points in fre- 
quency and momentum space. One can take as indepen- 
dent variables in T(k, w^) the modulus of k and positive 
frequencies uik, and we have adopted accordingly a grid 
of dimension I x I covering those variables. The advan- 
tage of this approach is that the number I plays the role 
of cutoff, making straightforward to assess the effect of 
higher frequencies and momenta as I is increased. We 
have solved the integral equation in I x / grids with I run- 
ning up to a value of 200, for which it is still viable to 
invert a matrix of dimension I 2 . One can rely moreover 
on the scale invariance of the theory to find the trend at 
larger values of I, as the critical coupling a c must obey 
a well-defined finite-size scaling law as a function of the 
cutoff a c (l) — a c (oo) + cjV . 

For a given value of N, we have determined the critical 
point at which the vertex r(k, uik) blows up. The value 
of a c (l) for our largest I provides an upper bound for the 
critical coupling in the continuum limit, as a c (l) turns 
out to be always a decreasing function of the variable 
I. On the other hand, the use of the above finite-size 
scaling law allows to estimate a c (oo) . We have cho- 
sen to represent in Fig. [2] the conservative upper bound 
a c (200) as a function of N. In marked difference with 
other approaches, we observe that now a critical cou- 
pling always exists, no matter how large the value of N 
may be. For N = 4 corresponding to graphene, we get 
the values a c (200) » 1.08 and a c (oo) w 0.99. 

We see that the value of the critical coupling obtained 
upon dynamical screening of the interaction is substan- 
tially smaller than the value found for N — 4 in the 
static approximation. This agrees with the results ob- 
tained in Ref. [2(| , where the resolution of the gap equa- 
tion was accomplished taking into account the frequency- 
dependent polarization. The critical coupling obtained 
there for N = 4, a c ~ 0.92, is actually very close to our 
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estimate a c (oo) ~ 0.99. These values are also close to 
the critical coupling a c « 1.08 found in lattice Monte 
Carlo simulations [14j. suggesting that the consideration 
of the dynamical screening provides a most sensible ap- 
proximation to the excitonic instability. 

Conclusion. — We have shown that the various many- 
body corrections used to dress the electron quasiparticles 
and the Coulomb interaction can have significant impact 
on the excitonic instability. The effects of the electron 
self-energy corrections and the electron-hole polarization 
have an important role in reducing the strength of the 
instability. We have seen however that, as anticipated in 
Ref. [20| . the simple static approximation overestimates 
the screening effects, and that a more accurate approach 
to the problem requires the consideration of the dynam- 
ical screening of the interaction. 

It is puzzling that, if we were to use the nominal values 
of the parameters appropriate for graphene, we would ar- 
rive at the conclusion that an isolated free-standing layer 
of the material (for which a « 2.2) should be in the phase 
of exciton condensation. This is at odds with the absence 
of any experimental observation of a gap in suspended 
graphene samples. Our many-body analysis shows that 



the only possible relevant effects that have been dismissed 
are those related to the scaling of the quasiparticle pa- 
rameters. In this respect, the growth of the Fermi veloc- 
ity at low energies [24j can have a deep impact to prevent 
the excitonic instability [2^, [27|. This effect, expressed 
by the scaling law (flU)) at e = 0, has been already ob- 
served in experimental samples of graphene at very low 
doping levels [28|. It is quite plausible that, at the low- 
energy scales where the gap could develop (about three 
orders of magnitude below the scale of the high-energy 
cutoff), the scaling of the Fermi velocity may have driven 
the effective coupling to such small values that the exci- 
tonic instability cannot then proceed. This can be one 
more consequence of the nontrivial scaling properties of 
the theory of Dirac fermions, implying that the electrons 
in graphene approach a noninteracting regime as they 
are observed at energies arbitrarily closer to the charge 
neutrality point. 
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